SPARTA (7 Jul 2014)
# flow around a sphere on a uniform grid
# particles are created initially and input at boundaries at stream velocity

variable            x index 10
variable            y index 10
variable            z index 10

variable            n equal 10*$x*$y*$z
variable            n equal 10*40*$y*$z
variable            n equal 10*40*50*$z
variable            n equal 10*40*50*50
variable            fnum equal 7.33e+21/v_n

seed	    	    12345
dimension   	    3
global              nrho 7.03E18
global              fnum ${fnum}
global              fnum 7.33e+15
timestep 	    1.E-5
global              gridcut 0.1
global              surfmax 300

boundary	    o o o
create_box          -5.0 5.0 -5.0 5.0 -5.0 5.0
Created orthogonal box = (-5 -5 -5) to (5 5 5)
create_grid         $x $y $z
create_grid         40 $y $z
create_grid         40 50 $z
create_grid         40 50 50
Created 100000 child grid cells
  parent cells = 1
  CPU time = 0.0160089 secs
  create/ghost percent = 83.4778 16.5222

balance_grid        rcb cell
Balance grid migrated 0 cells
  CPU time = 0.215703 secs
  reassign/sort/migrate/ghost percent = 3.29581 0.406975 0.790959 95.5063

species		    ar.species Ar

mixture             air Ar frac 1.0
mixture 	    air group species
mixture             air Ar vstream 2500.0 0.0 0.0 temp 300.00

fix		    in inflow air all

collide		    vss air ar.vss

read_surf	    1 data.sphere
  602 points
  1200 triangles
  -1 1 xlo xhi
  -1 1 ylo yhi
  -1 1 zlo zhi
  0.100631 min triangle edge length
  0.00439601 min triangle area
  400 = cells with surfs
  4040 = total surfs in all grid cells
  22 = max surfs in one grid cell
  0.402525 = min surf-size/cell-size ratio
  99360 240 400 = cells outside/inside/overlapping surfs
  400 = surf cells with 1,2,etc splits
  995.862 995.862 = cell-wise and global flow volume
  CPU time = 0.158532 secs
  read/surf2grid/error/ghost/inout percent = 3.06558 10.6193 5.38823 67.8279 13.099
surf_collide	    1 diffuse 300.0 1.0
surf_modify         collide 1 1

create_particles    air n 0
Created 955103 particles
  CPU time = 0.184331 secs

fix                 3 balance 100 1.1 rcb part
collide_modify      vremax 100 yes

#compute             1 grid all ndensity mass u v w
#fix                 1 ave/grid 10 100 1000 c_1 ave running
#dump                1 grid 1000 tmp_flow.* id xc yc zc f_1

#compute             2 surf all press ke
#fix		    2 ave/surf 10 100 1000 c_2 ave running
#dump                2 surf 1000 tmp_surf.* id f_2

compute             g grid all n
compute             max reduce max c_g[1]

stats_style	    step cpu np nattempt ncoll nscoll c_max
stats		    100

# equilibrate flow for 1000 steps
# then benchmark for 1000 steps

run 		    1000
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 98.877 98.877 98.877
  grid      (ave,min,max) = 17.7766 17.7766 17.7766
  surf      (ave,min,max) = 0.0641327 0.0641327 0.0641327
  total     (ave,min,max) = 117.481 117.481 117.481
Step CPU Np Natt Ncoll Nscoll max 
       0            0   955103        0        0        0           11 
     100     6.328711   953815      316      282      118           67 
     200    14.142211   955271      399      356      124           71 
     300    22.289169   961413      452      406      128           68 
     400    30.656419   976316      488      438      118           83 
     500    39.085101   991300      479      418      126           73 
     600    47.567008   996135      490      437      118           71 
     700    56.088207   999249      500      438      113           73 
     800    64.632075  1000212      484      435      116           71 
     900    73.177618  1000649      481      429      137           69 
    1000    81.728391  1000977      471      424       96           78 
Loop time of 81.7284 on 1 procs for 1000 steps with 1000977 particles

Particle moves    = 984026644 (984M)
Cells touched     = 1099259304 (1.1B)
Particle comms    = 0 (0K)
Boundary collides = 0 (0K)
Boundary exits    = 2734942 (2.73M)
SurfColl checks   = 63687584 (63.7M)
SurfColl occurs   = 113255 (0.113M)
Collide attempts  = 14431483 (14.4M)
Collide occurs    = 11064796 (11.1M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 1.20402e+07
Particle-moves/step: 984027
Cell-touches/particle/step: 1.1171
Particle comm iterations/step: 1
Particle fraction communicated: 0
Particle fraction colliding with boundary: 0
Particle fraction exiting boundary: 0.00277934
Surface-checks/particle/step: 0.0647214
Surface-collisions/particle/step: 0.000115093
Collision-attempts/particle/step: 0.0146657
Collisions/particle/step: 0.0112444

Move  time (%) = 38.7673 (47.4343)
Coll  time (%) = 16.8096 (20.5676)
Sort  time (%) = 23.4415 (28.6822)
Comm  time (%) = 0.331122 (0.405149)
Outpt time (%) = 0.240594 (0.294382)
Other time (%) = 2.13829 (2.61633)

Particles: 1.00098e+06 ave 1.00098e+06 max 1.00098e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Cells:      100000 ave 100000 max 100000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
GhostCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
collide_modify      vremax 100 no
run 		    1000
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 104.027 104.027 104.027
  grid      (ave,min,max) = 17.7766 17.7766 17.7766
  surf      (ave,min,max) = 0.0641327 0.0641327 0.0641327
  total     (ave,min,max) = 122.63 122.63 122.63
Step CPU Np Natt Ncoll Nscoll max 
    1000            0  1000977      471      424        0           78 
    1100    8.5531101  1000121      482      433      139           75 
    1200    17.107294  1000755      465      430      115           74 
    1300    25.661691  1000509      481      427      107           75 
    1400    34.211475  1000412      498      434      116           78 
    1500    42.761815  1000205      511      451      124           80 
    1600     51.32178  1000609      496      440      117           71 
    1700    59.872093  1000448      484      429      101           67 
    1800    68.424453  1000647      475      425      106           69 
    1900    76.979192  1000434      489      434      108           74 
    2000    85.527936   999920      496      433      103           66 
Loop time of 85.528 on 1 procs for 1000 steps with 999920 particles

Particle moves    = 1003313778 (1B)
Cells touched     = 1121026695 (1.12B)
Particle comms    = 0 (0K)
Boundary collides = 0 (0K)
Boundary exits    = 2780944 (2.78M)
SurfColl checks   = 67995943 (68M)
SurfColl occurs   = 117374 (0.117M)
Collide attempts  = 15411978 (15.4M)
Collide occurs    = 11828431 (11.8M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 1.17308e+07
Particle-moves/step: 1.00331e+06
Cell-touches/particle/step: 1.11732
Particle comm iterations/step: 1
Particle fraction communicated: 0
Particle fraction colliding with boundary: 0
Particle fraction exiting boundary: 0.00277176
Surface-checks/particle/step: 0.0677714
Surface-collisions/particle/step: 0.000116986
Collision-attempts/particle/step: 0.0153611
Collisions/particle/step: 0.0117894

Move  time (%) = 40.5147 (47.3702)
Coll  time (%) = 18.1258 (21.1928)
Sort  time (%) = 24.1688 (28.2584)
Comm  time (%) = 0.336845 (0.393842)
Outpt time (%) = 0.246379 (0.288069)
Other time (%) = 2.13539 (2.49671)

Particles: 999920 ave 999920 max 999920 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Cells:      100000 ave 100000 max 100000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
GhostCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
